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^P . The decays B^ — > D*'^^ and B^ -^ D*^j can be used for an extraction of 

\Vub\- When the b and c quarks are nearly degenerate the rate for these modes 

—f , can be determined in terms of other observed rates, namely BB mixing and 

Q^' D* — > D'y decay. To this end we introduce a novel application of heavy quark 

Q^' and flavor symmetries. Although somewhat unrealistic, this limit provides us 

D . with a first estimate of these rates. 
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I. INTRODUCTION 

The extraction of the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix element Vub, 
theoretically and experimentally, stands out as one of the prominent challenges of particle 
physics. Both its magnitude (through decay rates) and its phase relative to other CKM 
elements (through CP violation) are of considerable interest. Nevertheless, its determination 
remains elusive, ultimately because the CKM matrix describes the mixing of quarks, whereas 
of course only hadrons are observed. In the case of Vcb, much of this impediment is overcome 
by the application of the heavy quark effective theory (HQET) [jl]J^, which in particular 
relates the strong interaction matrix elements of systems with b and c quarks, eliminating 
much of the strong interaction uncertainty. However, the u quark is by no means heavy, and 
it is notoriously difficult to separate its CKM and strong interaction couplings. 

This is true even for the experimentally clean semileptonic modes, B — > iriu and B -^ 
piu, although dispersive bounds on the shapes of the form factors help to restrict such strong 
interaction uncertainties 0. Inclusive semileptonic rates for B -^ Xu^u are theoretically 
under better control but are plagued, since | Kb] ^ iKibI, by the preponderance of 5 — i> Xcdp 
everywhere in kinematic space except in the endpoint region of maximal lepton energy, 
where a c quark cannot kinematically be produced; however, in this small region it has 
proved necessary to include hadronic model dependence 0. A technique involving invariant 
hadronic mass spectra for invariant mass below mo is under better theoretical control but 
requires neutrino reconstruction |Q. 

Hadronic modes, containing everywhere strong interaction uncertainties, are even more 
problematic. Even if these difficulties could be tamed, one would still encounter the mixing 
of different weak topologies. For example, B -^ ti^K~ may proceed through either h — >■ 
usu (oc V^feV^*s) or a penguin b —>■ sg —^ suu (oc VJ^^ts); the presence of spectator quarks 
complicates the analysis. 

One may probe Vub using the decays 5+ -^ -D*"'"7 (bu -^ cs'-y) and B^ — > D*^'^ (bu -^ 
cdj) P, collectively B~^ — >• DJ^'j. Even though the second decay is Cabibbo-suppressed 
compared to the first, the lower reconstruction efficiency associated with D*~^ compared 
to D*^ makes both processes worth studying. The hard, monochromatic photon in these 
decays, 2.22 and 2.26 GeV, respectively, provides a distinctive experimental signature. 

As shown in Fig. 1, the invariant amplitude for these decays at 0{G];') consists of only 
one weak topology, because four flavor quantum numbers change {AB = —1, AC = +1, 
A{S,D) = +1, AU = —1), thus flxing the quark couplings to the W boson.Q The photon 
may couple to any of the quark lines, although one expects the dominant contributions from 
radiation emitted by the light u or s. The diagram with the photon emitted from the W 
is suppressed, of O^G'jp). On the other hand, since gluons may pair-produce quarks, these 
decays may be sensitive to multiparticle intermediate states, for example B~^ —>■ D^K^ — »• 
D*s^l 101 ' ^^^ ^^^ weak topology remains unchanged. 

In this paper we study a theoretical limit in which the decays B^ —>■ DJ^'-f are calculable 
from flrst principles in terms of other measured quantities. We assume both the b and c 



^At 0{G'p), the nontrivial diagrams include the di-penguin b -^ sg, ug ^ c and a box diagram 
with the internal quarks crossed. 



quarks are heavy and nearly degenerate, Aqcd ^ ^b, Aqcd <^ ^c and Sm = nib — rric '^ 
Aqcd- We do not expect this hmit to be a good approximation to reahty. Our calculation 
is a starting point for further investigations on controlled approximations of this rate. 



FIGURES 




FIG. 1. The unique weak diagram topology at 0{Gp) for B'^ — > Dl^j. The diagram for 
B^ -^ D*^'y simply changes s to d. 

II. FORMALISM 

The dominant contributions to the amplitude for B^ -^ -^Ts)^ ^^^ ^^^ ^^'^S distance 
processes B^ -^ Dt,^ -^ DI'X'^ and B^ -^ B*^^ -^ ^*(s)'y- "^^^ premise of this calculation 
is that the matrix elements of B^ -^ -^(s)' -^(s) ^ -^M^' ^"^^ ^^ ~^ B*~^'~f can be related 
by heavy quark and SU{3) flavor symmetries to the rates for BB mixing and i5^^ — > D{s)l- 
Alternately, we expect that lattice calculations similar to those that determine the BB 
mixing matrix element ^j can be employed to compute that of B^ -^ Dt-. directly. 

For the initial calculation presented here, we work in the generally-low (GL) velocity 
limit 

6m = rrib - rric ^ Aqcd -^ mb, (2.1) 

which means that the four-velocity of each heavy particle in the process is assumed to 
remain constant to lowest order. This limit is more restrictive than the slow-velocity (SV) 
limit considered by Shifman and Voloshin |^, Aqcd "^ ^^^ "C m^. To indicate the quality 
of this assumption, we point out that the value of 7d* = vb ■ vd* is 1.45 for the strange 

case and 1.50 for the nonstrange case. On the other hand, letting uib = mb + A — 3A^/4mb, 
niD* = nic + A + X^/Arric, one finds that •jd* is parametrically small, 

m| + ml* 

Heavy quark symmetries may be used to relate the matrix element for B -^ D(^s) to the 
one for B ^ B only if the velocity of the D(^s) in the B rest frame is parametrically small, 
say 0{A/mb) or O {5m, /nib). It is therefore sufficient to assume the weaker SV limit for this 
part of our argument. However, in order to relate the photon emission processes to the 
real decays B* — > i?7 and D* -^ D'j, it is necessary to impose the GL condition to ensure 
softness of the photon, as we now argue. 

The matrix element of the electromagnetic current j^ between vector {D* or B*) and 
pseudoscalar [D or B) states is characterized by a single Lorentz invariant form factor 



{p\e\3,mp) = -tg{q')e^.uX.e*yY 



(2.3) 



where q = {p — p')- The form factor at g^ = governs the on-shell decay of the vector to 
the pseudoscalar, which is characterized by a transition magnetic moment, g{0) = efi [see 



below, Eq. (3J^)]. However, imphcit in this definition is that the mesons are on-shell states; 
if the transition occurs with one of the states being virtual, as is the case here, then the g^ 
argument of the form factor must be computed accordingly. The relevant scale over which a 
hadronic form factor changes is Aqcd- For example, in the non-relativistic potential quark 
model of mesons, the scale that dictates the behavior of the form factor g is the constituent 
quark mass ~ Aqcd- In B^ -^ ^*{t)l ^^^ form factor appears with argument q^ ~ —5m?, 
and in the GL limit one may approximate g{—5m?) ^ g{Q) = efi. 

The formalism of HQET combined with chiral symmetry, relevant to the physics of 
mesons containing heavy quarks, was developed in Refs. 0. It was extended to radiative 
processes, such as D* —>■ D'-f, in Ref. [l^, while the Lagrangian for heavy quark-heavy 
antiquark mixing appeared in [|11|] . We use ingredients from all of these works, but for brevity 
include only notation relevant to the present process. While convenient, the formalism is 
not really a necessary framework for our calculation. In fact, light mesons enter diagrams 
only at higher order in the chiral expansion. By stating our calculation in the language of 
this effective Lagrangian we are setting the stage for further investigations of, for example, 
the effects of the finite strange quark mass. 

Since chiral symmetry is one of the ingredients of the Lagrangian, we begin by including 
the octet of pseudo-Nambu-Goldstone bosons in the usual nonlinear form ^: 



e = exp(ffl//). 



(2.4) 



where 
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with / ~ 131 MeV. Under the chiral SU(3)l x SU(3)r transformation {L,R), 



(2.6) 



with U implicitly defined by the equality of these two forms. 

The ground state vector P\^ and pseudoscalar P^^ fields that destroy mesons of heavy 
quark fiavor A and light antiquark fiavor a are incorporated into the 4x4 bispinor 
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(^Z7m - PAal5) , 
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(2.7) 



Under SU(4) heavy quark spin-fiavor symmetry transformations 5* and SU(3)l x SU(3)r 
transformations [/, Ha transforms as 



HAa ^-^ Sa Hsb U a, 

H^''^U\h'''s^b^. (2.8) 

Defining the chiral covariant derivative and axial current by 

Kb = I {('d^i - id^i^)^, = -jd^^ab + 0{U'), (2.9) 

one finds that DC,, DC,^ , DH, and DH have precisely the same transformation properties 
under chiral symmetry as their underlying fields, and 

{A'^)\^UMA^ydU^\. (2.10) 

Thus, one obtains the heavy meson Lagrangian invariant under chiral and heavy quark 
symmetry with the minimum number of derivatives, 

C = -TtH^\v ■ D\HAb + gTTH^''HAb4\l5- (2.11) 

It should be noted that the kinetic term is canonically normalized when the field H has 2, 
rather than the usual 2M, particles per unit volume. 

Since lowest-order HQET integrates out heavy antiparticle degrees of freedom, it is nec- 
essary to include such fields explicitly when they can appear in the asymptotic states. One 



defines [^ 



jjAa ^ (^p^*Aa^M _ pAa^^-^ 0_ 



ff^^ = 7°if t^V- (2.12) 

For example, while the field P^a destroys a vector meson of flavor content Ad, the field 
p*Aa (jgg^j-oys one of flavor content Aa; the two are related by a chosen charge conjugation 



convention WB. One demands that Hag and H^"^, or if" and H^ , have the same trans- 



formation properties under heavy quark symmetry and chiral transformations. Then the 
construction of the Lagrangian for mesons with heavy antifiavors is straightforward. 

As for direct mixing between the two sectors, the standard model operator that destroys 
quarks of flavors A, a and creates quarks of flavors B, b, given charge constraints —Qa + Qu = 
Qb - Qb, is 

C,ab,AB ^ ^^M(i _ ^^)a . 5^^(i _ ^^)5. (2.13) 

For BB mixing, one sees that A and B are b quarks, and a and b are d quarks; in the present 
case. A, B, a, and b are b, c, u, and [d or s) quarks, respectively. In the limit that HQET 
and SU(3) are good symmetries, both processes have the same strong matrix element. In 
Ref. |ll[| it was shown that ( |2.13|) matches in the symmetry limit of the effective theory onto 
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(2.14) 



where C is the color vacuum saturation coefficient, 1 for BB mixing and 3/8 for BD(s)- 

Calculation of the process also requires understanding the coupling of heavy mesons to 
photons. Minimal substitution for the kinetic term in Eq. ( p.ll|) provides no coupling to 
real transverse photons at lowest order, while minimal substitution applied to the second 
term in ( |2.11| ) requires including at least one pion for a nonvanishing term. The unique 
lowest-order electromagnetic term |]TU| is the transition magnetic moment operator 



5C 



-e/i^^Tr H HAaO''^'^^™ + charge conjugate. 



(2.15) 



III. CALCULATION 



The photon coupling in Eq. (|2.15| ) represents the lowest-order heavy-to-heavy electro- 
magnetic transition. It does not couple pseudoscalars to each other, but provides for the 
decay of vector mesons: 



r(p 



Aa 



PAal) = -a\^Aa\^E^. 



(3.1) 



To lowest order in HQET and SU(3), one may writeQ //^^ = Qa/^- For P\^ = D*^, the 
current experimental bound |l^ of r{D*^) < 0.131 MeV and electromagnetic branching 
fraction < 3.2%, lead to |/i| < 2.5 GeV^^ 

The calculation of the rate is then straightforward in the combined HQET and chiral 
symmetry limit, and consists of the two diagrams depicted in Fig. 2. We use the GL limit 
to compute the amplitude, but retain full mass dependence in the phase space calculation. 
The result is 
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(3.2) 



B 



with Vcs — ^ Vcd and m^,^ -^ m*j^ for B+ -^ D*+'j. 





FIG. 2. The two leading diagrams for B~^ — > D*,\^ that occur in the heavy quark and chiral 
limits. The double line indicates a heavy meson, one containing a b (c) quark on the left (right) of 
the flavor-changing vertex, indicated by a box. 



^In Eq. (14) of Ref. 0, M is labeled (3 



The effective Lagrangian mixing parameter (3 in ( |2.14|) is related to familiar quantities 
of BB mixing via 
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-jlniBBB = ^13. 



(3.3) 



A full treatment requires including renormalization point dependence in Bb-, but we neglect 
this here for simplicity. Now there are two obvious directions of analysis. One is to eliminate 
P using its calculated value from lattice simulations. Using |Q /^ = 170 ± 35 MeV and 
Bb = 0.98 ± 0.06 (renormalization point 2 GeV), and tub = 5.279 GeV, r^ = 1.6 ■ lO'^^ g^ 
one finds 4/3 = 0.40 ± 0.17 GeV^, and 
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(3.4) 



with Vcs -^ Vcd for 5+ -^ -0*^7, in which case the coefficient (taking \Vcd\ = 0.22) is 7 ■ 10~^. 
The other direction for analysis eliminates /? as a systematic strong interaction uncertainty 
between this process and BB mixing. In particular, Atub oc Gl^lV^f^V^^^p/?, which means|^ 
that eliminating P leads to values for the ratio |V^JbV^sP/l^6^dl'') which is proportional to 
(p^+?7^)/[(l— p)^+77^]^ in the usual Wolfenstein space. Curves of this family tend to intersect 
the p axis with vertical slope, which is useful since the current experimentally allowed region 
in p-i] space is broad in the p direction. 

IV. PROSPECTS 

The inclusion of HQET-violating, chiral symmetry-violating, and short-distance correc- 
tions is straightforward, and this program should certainly be carried out. Indeed, many of 



the necessary corrections already appear in the literature, including corrections to pAa |]T0 
meson decay constants 
corrections [p],! 



Tl| , p!3| , BB mixing [|TI|, and a number of other SU(3)-symmetry 
However, a number of other diagrams remain to be computed, because 
the inclusion of the photon changes the allowed quantum numbers of intermediate states. 
For example, there are new chiral loop diagrams for BD(^s^ mixing compared to BB mixing 
where a single pion emerges directly from the mixing vertex. 

More urgent and less tractable is understanding of the radiative form factor glq^) at re- 
alistic momenta. Taking g{q'^) — >• g{0) = ep introduces large uncertainties into our estimate 
for the decay rate. Alternatively, one may use only the less extreme SV limit to relate the 
mixing amplitudes and express the rate in terms of the corresponding form factors. This 



is done by replacing ep in Eq. (|3.2| ) by gBiq ) + gnil ) at q = —6m, which can be used 
a starting point for a calculation based on estimates using hadronic models for the form 
factors. 



^Note that T{B+ 



DJ^^) oc Gpf^"^, while Atub oc Gpf], due to the different weak topologies. 



We expect that the best determination of Vub will still require careful analysis of a number 
different decay modes; however, the decays B^ -^ -^M^ provide an important additional 
handle on the problem. 
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